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We construct a new order parameter for finite temperature QCD by considering the quark conden- 
sate for U(l)-valued temporal boundary conditions for the fermions. Fourier transformation with 
respect to the boundary condition defines the dual condensate. This quantity corresponds to an 
equivalence class of Polyakov loops, thereby being an order parameter for the center symmetry. We 
explore the duality relation between the quark condensate and these dressed Polyakov loops numer- 
ically, using quenched lattice QCD configurations below and above the QCD phase transition. It is 
demonstrated that the Dirac spectrum responds differently to changing the boundary condition, in 
a manner that reproduces the expected Polyakov loop pattern. We find the dressed Polyakov loops 
to be dominated by the lowest Dirac modes, in contrast to thin Polyakov loops investigated earlier. 
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Introductory remarks 

Understanding the nature of confinement has been a chal- 
lenging task for many years. Several scenarios with differ- 
ent candidates for the relevant gluonic excitations were 
proposed, but no closed picture has emerged yet (it is 
even still debated whether confinement is predominantly 
an infrared or an ultraviolet phenomenon). Also a con- 
nection of confinement to chiral symmetry and its break- 
ing has been conjectured, but not been shown either. 

In recent work [J [2i 9 have explored the idea of 
connecting quantities sensitive to confinement to spectral 
sums for Dirac and covariant Laplace operators. These 
ideas were developed further in 0, [^, Q where it was 
shown that also other quantities such as quark propaga- 
tors and heat kernels may be turned into order parame- 
ters for the breaking of center symmetry. Spectral sums 
provide a natural decomposition into infrared (IR) and 
ultraviolet (UV) parts and allow one to analyze their re- 
spective role in confinement, as studied numerically using 
quenched 0, d, IE £^nd dynamical 0] lattice config- 
urations. 

In this letter we build on those results and develop a 
new order parameter for center symmetry. In particular 
we Fourier transform the quark condensate (that turns 
into the chiral condensate in the massless limit) with 
respect to a U(l)-valued temporal boundary condition 
for the fermions, which we parameterize with a phase 
(fi £ [0, 27r). We show that this duality transformation 
turns the quark condensate into the expectation value of 
an equivalence class of Polyakov loops which all have the 
same winding number n e The winding number n is 
the conjugate variable to the phase ip. To the equivalence 
class of loops with winding number n = I, which trans- 
forms under center transformations in the same way as 
the conventional thin Polyakov loop, we refer to as the 
" dressed Polyakov loop" . 

Since for pure gauge theory the deconfinement transi- 
tion can be understood as spontaneous breaking of the 
center symmetry Q , the dressed Polyakov loop is an or- 
der parameter for confinement in pure gauge theory. The 



center transformation property of the dressed Polyakov 
loop is independent of the quark mass parameter m. It is 
particularly interesting that the mass parameter m may 
be used to relate the chiral condensate and the conven- 
tional thin Polyakov loop. In the limit where the quark 
mass parameter m is sent to infinity, the dressed Polyakov 
loop reduces to the thin loop that winds once. Con- 
versely, for the more interesting limit m — > we recover 
the dual of the conventional chiral condensate. We stress 
that this dual has a vanishing value below Tc, while it 
acquires a non-zero value above the transition. 

The dressed Polyakov loop can be written as a spectral 
sum for Dirac eigenvalues. Using quenched SU(3) gauge 
configurations below and above the critical temperature, 
we show that dressed loops are IR dominated (in con- 
trast to conventional thin Polyakov loops 0) and that 
the deconfinement transition is signaled by a changing 
dependence of the IR modes on the boundary condition. 

Dual quark condensate 

We begin our discussion with recalling the definition of 
the chiral condensate E. The starting point is the scalar 
expectation value S(to, V) of the fermion bilinear ipip 
evaluated at finite volume V and mass m, 



S(m,y) = - 



V 



(1) 

In the second step the fermions were integrated out and 
the remaining expectation value {■■■)g is the path inte- 
gral over the gauge fields with gauge action and fermion 
determinant included in the weight factor. D denotes the 
Dirac operator at vanishing quark mass. We refer to ([T]) 
as the "quark condensate". The chiral condensate, i.e., 
the proper order parameter for chiral symmetry break- 
ing is obtained through a double limit, where first the 
4-volume V is sent to infinity and then the quark mass 
m to zero: S = lim,„_^o liniv-^oo S(m, V). 

An important result, that we will return to later, is the 
Banks-Casher formula @ which relates the chiral con- 
densate to the density p{0) of eigenvalues at the origin: 
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E = 7rp(0). Below the critical temperature Tc, where chi- 
ral symmetry is broken (E ^ 0), the eigenvalue density 
p(0) at the origin is non- vanishing, while above Tc the 
Dirac spectrum develops a gap, p(0) vanishes and chiral 
symmetry is restored (S = 0). Chiral symmetry break- 
ing is a feature of full QCD, i.e., massless fermions are 
taken into account. However, the above discussed mech- 
anism of a non-vanishing spectral density p(0) below 
and a spectral gap above Tc is known to hold also for the 
quenched case, i.e., pure gauge theory (see, e.g., |10|] for 
a lattice study of this property). 

We now define the dual quark condensate which we 
will later identify as an order parameter for center sym- 
metry. We work in a finite Euclidean volume with tem- 
poral extent (3 = l/ksT. For the fermion fields ip 
we use the generalized temporal boundary conditions 
4>{x,l3) = e"^ ■0(2^1 0). The canonical choice is anti- 
periodic, i.e., if = TT, while here we allow for arbitrary 
values (fi G [0, 27r]. For the spatial directions the fermions 
obey periodic boundary conditions, and the gauge fields 
are periodic in all four directions. 

Now the quark condensate is considered for an arbi- 
trary boundary angle ip indicated by a subscript for the 
Dirac operator. The "dual quark condensate" S„ is de- 
fined as the Fourier transform with respect to tp, 



I]„(m, V) 



27T 



J9:(t.[(™ + I>J-])^, (2, 



where the index n is an integer. This gauge invariant 
quantity, in particular the case n = I, will be investigated 
analytically and numerically now. 

Relation to dressed Polyakov loops 

Using the lattice rcgularization, we now discuss the re- 
lation of the dual condensate to equivalence classes of 
Polyakov loops (dressed Polyakov loops). To be specific, 
we use the staggered lattice Dirac operator (other lattice 
Dirac operators give equivalent results) 



2a 



(3) 



with the staggered sign function ii^{x) = (— 1)^"^+ ■■■ . 
The coordinates x, y run over all sites of a 4-dimensional 
X N4 lattice with lattice spacing a. The gauge link 
variables Ufj,{x) are elements of the gauge group SU(N). 

The U(l)-valued temporal fermionic boundary condi- 
tions are most conveniently introduced by attaching the 
boundary phase to the temporal link on the last time- 
slice, U4{x, N4) e**^ U4{x, N4) (at the moment ip is held 
fixed and the Fourier transformation of ^ is performed 
only later in Eq. Inserting this into the Dirac op- 

erator, we evaluate the propagator in ([2]) for sufficiently 
large m as a geometric series 



Tr[im + D,n . ^ £ ^Tr [(D,)'=] 



fe=o 



(4) 



The Dirac operator ([3]) contains only terms that connect 
nearest neighbors. The power (Di^)*^ corresponds to a 
chain of k hops. The trace in (U is over color- and space- 
time indices. The latter trace implies that the chains of 
hops have to form closed loops I (thus on a lattice with 
even numbers of sites in all directions k must be even). 
Consequently, the sum in ^ can be reexpressed as a sum 
over the set C of all possible closed loops on the lattice. 



Tr[(m + D^)-i] 



1 



s{l)e 



lipqil) 



m ^ (2am)l'l 



-^^c^U^ix). (5) 



The remaining trace Trc is over the color indices of the 
ordered product of all link variables Uf_t{x) in a loop /, and 
we use U-fj,{x) = U^{x — fiy . By s{l) we denote the sign 
of a particular loop I which is obtained as product of the 
staggered sign factors. Each step in the loop comes with 
a factor of l/2am from the discretization in ^ and the 
normalization in The number of steps, i.e., the length 
of the loop is denoted by |/|. The loops may close around 
the boundary. When they close around the temporal 
boundary, they pick up a factor of exp(iip) if they run 
forward in time, and a factor of exp{—iip) for backward 
running. Denoting the number of times a loop I winds 
around the compact time direction by its winding number 
q{l) e we obtain the factor exp{i(pq{l)) in ([5]). 

This is how the boundary angle ip can be used to dis- 
tinguish between closed loops of different winding number 
0]. When expression (O is inserted into the formula ^ 
for the dual condensate, the 93-integration with the ad- 
ditional Fourier factor exp(— j(^?i) projects to loops of a 
particular winding number n. We finally obtain 



;e£(n) ^ ' 



Tr, n U,ix)\ , (6) 



where the sum now runs over the set C{n) of loops that 
wind n-times around the compact time direction. The 
case of n = 1, i.e., the dual condensate Y,i{m^V) which 
corresponds to loops that wind exactly once, is what we 
refer to as the "dressed Polyakov loop". From ^ it is 
obvious that in the large-m limit the dominant contribu- 
tion is the conventional thin Polyakov loop (as this is the 
shortest loop winding once). 

Let us finally discuss the behavior of the dual con- 
densates under a center transformation, {74(0;, to) ~* 
zU4{x,to), where all temporal links on a time-slice, i.e., 
at some fixed to, are multiplied with an element z of the 
center of the gauge group. The gauge action and the 
measure are invariant under this transformation, but the 
center symmetry is broken spontaneously at the critical 
temperature Tc, signaling the transition to the deconfined 
phase ^8]. Both the thin Polyakov loop and the dressed 
Polyakov loop are order parameters for center symme- 
try. More generally a loop that winds n-times around 
the compact time direction picks up a net factor of 
and thus our dual condensajtes transform under a center 
transformation as I]„ — > 2;"E„. In particular the dressed 
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FIG. 1: The integrand 7((^) = F-iE,(('" + A<:')-')g of® 
in lattice units for two values of am. The data are from 20 
gauge configurations on 12'^ x 6 lattices below (T — 255 MeV, 
a = 0.129 fm) and above Tc (T = 337 MeV, a = 0.098 fm). 

Polyakov loop, i.e., the dual condensate for n = 1, trans- 
forms as El — > zT,i, which is the same transformation 
law as for the thin Polyakov loop. 

Spectral sums and numerical analysis 

An elegant way to express the dual condensate ([2]) is 

as a spectral sum over all Dirac eigenvalues A^''* (again 
evaluated for boundary angle ip), 

i 

First of all this spectral representation consists of a finite 
sum over the eigenvalues (due to the regularization on a 
finite lattice) , whereas the representation ([6]) contains an 
infinite sum over loops (even on a finite lattice). 

Secondly, the sum on the rhs. of ([7]) can easily be com- 
puted numerically and then allows for the appealing pos- 
sibility to study the role of individual parts of the spec- 
trum. Since the (purely imaginary) eigenvalues appear 
in the denominator, we expect the dual condensate to be 
dominated by IR modes (for small mass m) and thus to 
have a well-defined continuum limit (see also the discus- 
sion of the convergence of the spectral sums in [^). We 
stress that the above discussed properties under center 
transformations are independent of the mass parameter 
TO. Naturally one is interested in small values of m, where 
the spectral sums are IR dominated. For the Figures 1 
and 2 below we use masses near the strange quark mass 
and only in Fig. 3, where we illustrate the transition to 
the large mass behavior, we increase the mass to 771 = 1 
GeV. 

For our numerical analysis of the spectral representa- 
tion ([7]) we compute complete spectra of the staggered 
Dirac operator ^ at different boundary angles (p with 
parallel LAPACK routines. We use quenched SU(3) 
gauge configurations generated from the Liischer-Weisz 



FIG. 2: The dressed Polyakov loop at m = 100 MeV in units 
of GeV"^ as a function of the temperature T in MeV. 

gauge action [l3|. The scale was determined [l^ from 
the Sommer parameter setting tq = 0.5 fm. We use 
X N4 lattices with L ranging from 8 to 14 and N4 
from 4 to 8, and adjust the couplings such that we have 
ensembles below and above the critical temperature Tc ~ 
300 MeV. For the gauge configurations above Tc the thin 
Polyakov loop has a non- vanishing expectation value and 
for our numerical study we use configurations where the 
thin Polyakov loop is real. Below we discuss how the 
spectral sum for the dressed Polyakov loop may generate 
also the two possible complex phases. All errors shown 
are statistical errors from the Jackknife method. 

We begin with analyzing the (/^-dependence of the in- 
tegrand in ([T]) . In Fig. [T] we plot this integrand versus 
for two values of am comparing an ensemble below Tc to 
one above Tc- It is obvious that below Tc the integrand 
is essentially constant, while above Tc it shows a pro- 
nounced cosine type of behavior. Integrating over ip with 
the weight exp(— i(^), i.e., n = I, gives a vanishing dressed 
Polyakov below Tc , while above Tc a non- vanishing value 
is observed. We conclude that the transition from con- 
finement to deconfinement leads to a different response 
of the spectral sums to the changing temporal fermion 
boundary conditions. Below we will demonstrate that 
the IR modes play the dominant role in this process. 

Above Tc, when using instead of an ensemble with 
real thin Polyakov loop, configurations where the thin 
Polyakov loop has one of the two complex phases, the in- 
tegrand is shifted by ±27r/3, but otherwise has the same 
form as depicted in Fig. [TJ Integrating this shifted in- 
tegrand shows that the dressed Polyakov loop produces 
the same phase pattern as the thin Polyakov loop. 

Now we demonstrate that the dressed Polyakov loop 
does indeed signal the phase transition. In Fig. [5] we show 
the results for the dressed Polyakov loop aX m — 100 
MeV as a function of T. The necessary (/^-integration 
was implemented with the extended Simpson rule using 
typically 8 or 16 values of f. It is obvious that below 
Tc the dressed Polyakov loop vanishes, while above Tc it 
assumes a non- vanishing value, signaling that the center 
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FIG. 3: Tlie upper plots show the (normalized) individual 
contributions jC(|A|)/Ei| to ((7|) versus |A| for two values of 
m, while the lower plots are for the normalized accumulated 
contributions |^(|A|)/Ei| as a function of |A|. The data are 
from the same ensembles already used in Fig. 1. 



procedure, we compare two different values of to. In both 
cases the largest contributions come from the IR end of 
the spectrum. However, for the smaller mass the contri- 
butions beyond the deep IR have died out completely. 

When considering the relative role of IR and UV 
contributions, one must take into account that the 
density of eigenvalues increases strongly with increasing 
|A|. Thus in the lower plots of Fig. [3] we show the 
accumulated contribution j4(|A|) — J2\x'\<\x\ ^i^')- 

particular we plot the normalized function |^(|A|)/S]i| 
which approaches 1 in the UV limit. Here the differ- 
ence between the large and the small mass is rather 
pronounced: For the large mass, where the dressed loop 
is closer to the thin Polyakov loop, the accumulated IR 
modes overshoot the value of 1 and the UV contributions 
are necessary to get the correct value. This is similar 
to the observations in 0, i, 0] for the thin Polyakov 
loops. Towards the chiral limit, however, the approach 
to the final value of 1 becomes monotonic, and only the 
IR modes give sizable contributions to the spectral sum 
for the dressed Polyakov loop. The scale up to which 
eigenvalues are relevant grows with T. 



symmetry is broken in the deconfining phase. 

Another interesting connection between quark conden- 
sate and dressed Polyakov loop is obtained by using the 
Banks-Casher type of representation already addressed 
above. After performing the consecutive limits of infinite 
volume and vanishing mass, the chiral condensate can be 
written as the density of eigenvalues at the origin also for 
arbitrary boundary angle if. The dual condensate is then 
obtained by integrating the i^j-dependent spectral density 
p{0)ip, and for the case of n = 1 we find 

Si = l'^^e-^^p{0)^. (8) 

Below Tc the spectral density at the origin is constant 
as a function of ip and a vanishing dressed Polyakov 
loop emerges. More interesting is the situation above 
Tc, where a non-trivial ip dependence is necessary for a 
non- vanishing dressed Polyakov loop. Naively one would 
think that above Tc the spectral density at the origin 
must be zero, such that the chiral condensate may van- 
ish. However, in [TTl (for a different phase convention) it 
was shown that the spectral gap depends on the relative 
phase between the boundary angle if and the phase 9 of 
the Polyakov loop. If (p equals the negative Polyakov loop 
phase the gap closes completely, giving rise to a non-zero 
spectral density. Inserting p{0)ip oc S{(p + 6) in ^ one 
obtains a non-vanishing dressed Polyakov loop above Tc 
with the correct phase 6. 

Let us now consider the individual contributions 
C(A) = i2nV)-^ J dip exp{-itp){{m + \^)-^)^ to the 
spectral sum ([7]). The upper plots of Fig. [3] show 
|C(A)/Ei| as a function of |A|. Since the case of vanishing 
quark mass to can be obtained only through a limiting 



Concluding remarks 

In this letter we have shown that a duality transforma- 
tion of the quark condensate with respect to the fermionic 
temporal boundary condition gives rise to an order pa- 
rameter for center symmetry. This order parameter can 
be viewed as a set of closed loops with the same wind- 
ing number around compact time. For the case of single 
winding, which we refer to as the dressed Polyakov loop, 
our observable interpolates between the quark conden- 
sate (via a Fourier transform) and the thin Polyakov loop, 
in the limits of vanishing and infinite mass, respectively. 

We studied the corresponding spectral sums of Dirac 
operator eigenvalues numerically for quenched gauge con- 
figurations. It was shown that the transition from the 
confined to the deconfined phase is seen as a different de- 
pendence of the spectral sums on the fermionic boundary 
condition. Decomposing the spectral sum shows that the 
main signal comes from the IR part of the spectrum. 

As an outlook we stress that our theoretical consid- 
erations also hold for other gauge groups (where, e.g., 
for SU(2) the phase transition can be of different order) 
and for the case of dynamical quarks, the consequences 
of which would be interesting to study numerically. 
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